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Abstract 

Physics beyond the standard model (SM) can be parameterized with an effective La- 
^^ ' grangian that respects the symmetries of the standard model and contains many operators 

of dimension six. We consider the subset of these operators that is responsible for flavor di- 
agonal anomalous color magnetic (CMDM) and electric (CEDM) dipole couplings between 
-^ quarks and gluons. Invariance of these operators under the SM implies that they contribute 

to Higgs boson production at the LHC and we study the corresponding constraints that 
can be placed on them. For the case of the top-quark we first review constraints from 
top-quark pair production and decay, and then compare them to what can be achieved by 
i-£h . studying tth production. We also constrain the corresponding couplings for 6-quarks and 

light quarks by studying pp — > bbh and pp — > hX respectively. 
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I. INTRODUCTION 

The flavor diagonal dipole couplings between quarks and gluons of magnetic 
and electric type are generalizations of the electric dipole moment and anomalous 
magnetic moment of fermions given by the dimension five operators, 

£ = | / TV" (a? + i 78 df) / G%,. (1) 

The operator in Eq. [I] is not gauge invariant under the full standard model gauge 
group and this needs to be remedied in the context of effective field theories which 
provide a general framework to study physics beyond the SM. Within this formalism 
one assumes that there is some new physics (NP) that appears at a high energy scale 
A that can be integrated out and represented by an effective Lagrangian valid for 
energies below A. The dominant effects of the new physics at scales below A are 
then described by the lowest dimension operators, and complete catalogs for these 
operators (at least up to dimension six) exist in the literature [l|, |2[. 

After last year's discovery of a light Higgs boson with mass rrih ~ 126 GeV 1 , the 
gauge invariance of the effective Lagrangian must be imposed within the framework 
of the linear realization of electroweak symmetry breaking. The effective Lagrangian 
thus includes the Higgs field as an active degree of freedom. This has implications for 
the operators in Eq. [1] they really are of dimension six and they contain couplings 
to Higgs bosons. The explicitly gauge invariant form is given by (in the notation of 
Ref. (J) 

C = 9s^§ qa^T a u ^ + g*fjjL qcr^T a d <PG% + h.c. (2) 

where q is the SM quark doublet, <fi is the scalar doublet, <pi = eij<pj and the 
SU(3) generators are normalized as Tr(T a TJ,) = 5 at b/2- Electroweak symmetry is 
spontaneously broken when the scalar acquires a vacuum expectation value (vev) 
<(/)>= v/y/2, v w 246 GeV. An expansion of this gauge invariant Lagrangian 
reproduces the terms in EqJT] with the identification 

9 V2v 

a u,d = —£T Re ( d uG,dG> 

d u,d = -jr M<Wg)- (3) 

The effect of imposing invariance under the SM group is then that the operators 
of Eq. [I] can only exist in conjunction with corresponding operators that involve a 
Higgs field. This opens up the possibility of constraining the couplings by studying 
processes with a Higgs boson, such as pp — > ffh which we do in this paper. In all 
cases, when we quote constraints on d q c we will implicitly assume a new physics 
scale A = 3 TeV. 



1 Interpreting this state to be an elementary scalar. 



II. CASE OF THE TOP-QUARK 

The flavor diagonal dipole couplings of the top-quark to gluons are known in 
the literature as the CEDM and the CMDM and have been studied at length in 
connection to top-quark pair production and decay. The notation, however is not 
uniform and different definitions exist |6l-ll5|. A typical result is that of Ref. [14j 



where it is found that using CP-odd observables the 5a statistical sensitivitv with 
10 lb-« to d' t is of order 0.1/m,. The CMDM has also been studied before fllB- 



18] and constraints have been derived from production cross-sections. Very recently 



Ref. [19] has found that the study of spin correlations over the life span of the LHC 
can produce limits at the level of 0.03/m 4 for both a? and df . 

We begin this section revisiting the constraints from top-quark pair production 
cross-sections and T-odd asymmetries both at 8 and 14 TeV. These will serve as 
comparison points for our study of constraints from tih processes, a possibility that 



has also been known for a long time 20| and has been revisited recently 23|, |24 



Our results for the top-quark couplings are in agreement with these recent papers 
and we expand on them by considering the T odd correlations that would allow one 
to separate the effect of a CEDM from that of a CMDM and by extending the study 
to include light quarks. 



A. tt production and decay 

We first consider ti production at the LHC at 8 TeV in order to compare with 
recent experimental results. Both CMS [3| and ATLAS have results that have not 
yet been combined. Since all the existing results are compatible, we will use one 
of them to illustrate the constraints on the new couplings. Taking the ATLAS 
cross-section in the lepton plus jets channel [5] with errors added in quadrature (but 
dominated by systematic error) 

a(tt) = (241 ± 32) pb, (4) 



and comparing with the theoretical expectation obtained from HATHOR [4j as 
quoted by ATLAS, 

a(tt) = (238+^) pb (5) 

we find a representative ratio between the measured and expected cross-sections 

ff{tt)Exp 



a{tt) 



1.01 ±0.17. (6) 



TH 



In this section we will interpret the error in this ratio as the room that remains for 
new physics that can affect top-quark pair production. 

We begin our calculation by implementing the Lagrangian of Eq. [2] into 



MadGraph5 25| with the aid of FeynRules 26|. We use the resulting model UFO 



file to generate top-quark pair events for different values of af and d 9 in a range that 
changes the SM cross-section by factors of a few. We then fit the numerical results 
to obtain approximate expressions for cross-sections and asymmetries in terms of 
the new couplings. The results of our simulations are presented in the Appendix in 
Figure El and the corresponding fit to these points is given in Eq. IA1I 

In Figure [I] we compare the fit to the current measurement in the form of Eq. [6l 
and from this comparison we extract the current LHC exclusion region. The central 
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FIG. 1: dtG parameter space allowed by the measurement of the it production cross-section 
at the LHC at 8 TeV. The center curve (blue) corresponds to the central value of Eq. [HJ 
and the shaded region is that allowed at la. 

value (for which the measured cross-section is approximately equal to the SM) is 
reproduced along the curve that goes through the origin. Along this curve there is 
a cancellation between pure NP contributions and interference between NP and the 
SM. The shaded region shows the parameters allowed by Eq. [6] at the la level. 

Constraining these couplings one at a time, results in the la ranges for the top 
CMDM and CEDM 



0.034 < m t a 9 t < 0.031 or -0.55 < m t a 9 t < -0.46, 

and \d 9 t | < 0.12/mt. 



(7) 



From Figure [TJ we see that there are two allowed ranges for the real part of d t G 
(proportional to af). Of course the effective Lagrangian formalism assumes that the 
new physics is small compared to the SM and this makes the range closest to the 
origin, presented in Table HI the more natural one. Approximately, this is the range 
one would get by keeping only the term linear in the CMDM in Eq. [6j This figure 
is in agreement with that of Ref. 13[ with an opposite sign convention for dtG- 

Next we consider possible improvements at 14 TeV. It should be clear that simply 
measuring the ti production cross-section will not change the above picture much 



as long as the uncertainty in the cross-section is dominated by systematics. More 
interesting is the possibility of using large statistics samples at 14 TeV to measure 
asymmetries in specific channels. In particular T-odd correlations are linear in the 
CEDM and potentially have a better sensitivity to this type of new physics. Here 
we will focus on the clean di-muon channel but generalizations to other channels 
are possible and have been discussed before 32j. For our numerical study we use 
default MadGraph5 cuts that require the intermediate top-quarks and W bosons in 
the process pp — > ti — > bbp + pTvv to be within 15 widths of their mass shell. We use 
MSTW2008 parton distribution functions and a set of acceptance and separation 
cuts for muons and b quarks given by 

Pt(^) > 20 GeV p T (b,b) > 25 GeV 
177(6,6,^)1 < 2.5 AR(bb) > 0.4. (8) 

In addition we use a missing energy requirement, $ T > 30 GeV, that is known to 
significantly reduce the background. 

As detailed in the Appendix, we estimate the dependence of the cross-section 
on d t G by generating event samples at different values of the CEDM and CMDM 
and extrapolating with a fit, Eq. IA2I This fit is shown in Figure [2} where the 
central (blue) curve corresponds to a cross-section matching the central value of 
the SM prediction. The shaded region shows the parameter space for which the 
NP contributions are within the 1 a errors in the theoretical SM cross-section at 
next-to- leading-order (NLO), o 
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884llo? pb 27J Once again we can summarize the 
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FIG. 2: dtc parameter space for which the predicted tt production cross-section at the 
LHC at 14 TeV is within the la theoretical error in the SM cross-section at NLO. 

results taking only one non-zero coupling at a time as 

- 0.029 < m t a g t < 0.024 or -0.480 < m t a\ < -0.478, 



and \d 9 t \ < 0.10/m t . 



(9) 



As expected, this measurement will not result in a significant improvement over the 
bound with 8 TeV data, Eq. 0, unless the experimental systematic errors and the 
theoretical errors in the SM can both be reduced to the few percent level. As before, 
we only include the more natural range for of in Table HI 

Next we consider a T-odd correlation that can single out the CP violating CEDM 
(a correlation that has been studied before). As opposed to CP even observables 
such as cross-sections, this asymmetry is linear in the CEDM and therefore more 
sensitive to it. To quantify our bounds we consider the integrated asymmetry in the 
lab frame defined by 

a = *(o.>o)-*(o.<Q), (10) 

associated with the correlation 

Ox = Pt-(p^+ xp». (11) 

We choose Eq. [TTJ as this is the T-odd correlation most sensitive to the CEDM 



out of the many possible for this process [14J. One correlation will suffice for our 
purpose of comparing the bounds that can be obtained from processes with and 
without Higgs bosons. Of course, a more realistic discussion of this correlation 
requires using momenta that can be reconstructed; discussion of backgrounds; and 
possible combinations of different channels. All these aspects have been discussed 
previously in the literature and we do not repeat that discussion. 

Notice that we have defined the integrated asymmetry in Eq. [10] using the SM 
cross-section in the denominator. This choice ensures that the asymmetry is linear 
in the CEDM and does not depend on the new physics indirectly, through the 
cross-section. Experimentally it is simpler to define a counting asymmetry which 
would correspond to Eq. [10] with the cross-section as a function of the CEDM in the 
denominator. Such counting asymmetry is only approximately equal to A for small 
values of the CEDM. We present our numerical results for the asymmetry for four 
different values of Im(d t G), as well as a linear fit to these results in Figure [3] The fit 
indicates that a measurement with 10 fb _1 would have a (lcr) statistical sensitivity 
of 

|d?| < 0.009/mj, (12) 

an order of magnitude better than what can be achieved from a measurement of 
the cross-section, Eq. [9] This sensitivity is comparable to what can be achieved for 
a\ from cross-section measurements, reflecting the fact that the CMDM contributes 
linearly to the cross-section. In principle it is also possible to obtain increased 
sensitivity to a\ using other observables such as spin correlations, but existing results 
in the literature suggest the improvement is only by factors of a few [19| . 

The best existing indirect constraint on the top CEDM arises from contributions 
to the neutron electric dipole moment (edm) [20j (rescaled to the current limit |2l| ) 



via the Weinberg operator 22[ giving \m t d 9 t \ ^ 2.4 x 10 -4 . Although the LHC will 
probably not reach this limit, direct searches are always complementary to indirect 
constraints which rely on multiple assumptions. 
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FIG. 3: Asymmetries in the process pp — > tt in the di-muon channel as calculated with 
MadGraph5 for different values of the anomalous coupling dtc with the cuts of Eq. [5] and 
the corresponding fit. 

B. tth production 



We turn our attention to the Higgs production process in association with a 
top-quark pair. We begin as usual with a study of the cross-section which we 
calculate numerically at leading-order (LO) with MadGraph5. Our numerical results 
for several values of the top CEDM and the top CMDM couplings as well as a fit 
to these Monte Carlo (MC) points are presented in the Appendix in Figure [10] and 
Eq. IB1I We normalize our results to the SM cross-section and use them to estimate 
the sensitivity that can be obtained at V5 = 14 TeV. 

In Figure H] we illustrate the region in the a\ — d\ parameter space that would 
be allowed by assuming that the central value of a future measurement agrees with 
the SM value (blue curve). The shaded region corresponds to that in which the 
new physics effects modify the SM cross-section by at most la using the theoretical 
errors in the NLO cross-section for the SM, a(pp — > tih)NLO = (611^^ ) fb 28-130 



For tth we only consider the 14 TeV LHC because the existing 8 TeV results for this 



mode are only sensitive to cross-sections about 5-6 times larger than the SM 31 



Assuming that the experimental error can be kept below the theoretical uncertainty 
at NLO a future measurement can reach a 1-a sensitivity to the ranges 



0.06 < m t a 9 t < 0.03 or 
and 



0.016 < m t a\ < 0.008, 
< 0.02/m 4 . 



(13) 



As before, we have only included the more natural range for a g t of the two allowed 
solutions in Table HI We see a significant improvement over what can be achieved 
with the top-quark pair process, Eq. [9] for the CEDM and a modest improvement 
for the CMDM. 

Next we can try to isolate the CP violating coupling Im(d t G) with the use of 
the T-odd asymmetries A12, associated with the correlation 0\ of Eq. [TTJbut for 
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FIG. 4: dtG parameter space as it would be constrained by a measurement of the it/i 
production cross-section at the LHC at 14 TeV. The central (blue) curve corresponds to 
cross-sections that are equal to the central value of the SM. The shaded region is the one 
in which the NP contributions are below the la theoretical error for the SM at NLO. 



the process pp — y tth — y bbhfi + p uu, and with a second correlation 



O, 



Pbeam ' (p^+ X p^- ) p be 



(P/- 



Pm- 



(14) 



In principle this process admits additional correlations involving the Higgs momenta. 
Numerically, however, we find all of them to be significantly smaller. For our study 
of the asymmetries we assume that the top quarks decay as in the SM and consider 
only the di-muon channel in order to construct them, keeping in mind that the 
statistics can be increased by using other channels. No attempt is made to include 
the Higgs decay or the complications involved with its reconstruction. The results 
of our MC simulation for a few values of the top CEDM as well as a linear fit to 
those points are shown in Figure O Using the central value of the NLO cross-section 
for pp — y tth at 14 TeV in the SM, 611 fb 30], we see that to measure percent level 
asymmetries in this channel at the la level would require about a thousand fb -1 . 
When these integrated luminosities have been collected the asymmetries in Figure [5] 
could reach a sensitivity to m t d 9 t at the 0.007 level. Obviously one can improve the 
statistics significantly by considering other top decay channels, but the number of 
tth events is small so this type of measurement will not be possible any time soon. 



III. CASE OF THE BOTTOM QUARK 

We now turn our attention to the case of the 6-quark anomalous CEDM and 
CMDM couplings. In this case the bbh process is possibly the only handle on these 
couplings at the LHC as the bb cross-section is completely dominated by QCD and 
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FIG. 5: T-odd asymmetries with the largest sensitivity to the top quark CEDM in the 
process tth. The points correspond to MadGraph5 simulations for four different values of 
the CEDM, and the line is a linear fit to these points. 



the contributions of the new couplings are negligible. The production of bb pairs 
in association with a Higgs boson has been studied before, with a NLO result for 
m H = 120 GeV of J33 



a(pp -> bbhX) 



SM 



(5.8 ±1.0) x 10 2 fb. 



(15) 



More recently this mode has been studied in connection to Higgs bosons in the 
MSSM where the bbh coupling is enhanced for large tan/3 34 1. 

As before, we generate MC samples for several values of the CEDM and CMDM 
of the 6-quark chosen so that they change the SM cross-section by factors of a 
few, and we then perform a fit to these values. We generate the bbh events with 
default MadGraph5 cuts, without h decay, and without backgrounds. These results 
are presented in the Appendix in Figure [Tl] and Eq. IC1I respectively. 

Our fit can then be used to estimate the statistical sensitivity to the b-quark 
CEDM and CMDM that can be achieved with a measurement of the bbh cross- 
section. In Figure |6] we illustrate the region that would be allowed by requiring 
the new physics to remain within ler of the NLO theoretical cross-section. The 
central curve (blue) corresponds to values of dbc which reproduce the central value 
of the SM cross-section. We can extract from the figure the la ranges with only one 
coupling at a time (using m& = 4.2 GeV) 



- 1.3 x 10" 4 < m b a 9 b < 2.4 x 10~ 4 \d 9 b \ < (1.7 x 10' 



4 )M. 



(16) 



Unlike the case of the top-quark, the interference between the b-quark CMDM and 
the SM is very small so that its effect on the cross-section is dominated by a quadratic 
term, Eq. IC1I For this reason there is only one range for the allowed a 9 b in Eq. [16j 

The best existing indirect constraint on the b-quark CEDM arises from contri- 
butions to the neutron edm 20] via the Weinberg operator 22]. Rescaling to the 
current limit on the neutron edm, we find \m b d 9 b \ < 2 x 10~ 8 . 
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FIG. 6: dbc parameter space allowed by a measurement of the bbh production cross-section 
at the LHC at 14 TeV. The shaded area shows the region in which the NP contributions 
are below the la theoretical error in the NLO SM result. 

IV. THE LIGHT QUARKS INCLUDING CHARM 

For LHC processes involving only the light quarks, the effects of anomalous a 9 
or d 9 couplings are buried in QCD background and Higgs production offers a unique 
window into the color dipole operators at high energy. 

The light quark CEDM and CMDM couplings contribute to Higgs production 
at the LHC dominantly through the parton level process qg — > qh and to a lesser 
extent qq — > hg. These processes are dominated by the charm-quark within the SM. 
There are several implications from this observation: 

• We can look for the effect of the new couplings in the LO Higgs production 
process pp — > hX. Our strategy to constrain the new couplings in this case, is 
to require the new contributions to fall within the theoretical uncertainty of 
the dominant SM gluon fusion cross-section and we pursue this route below. 

• The SM tree-level contribution to Higgs production from qg — y qh and qq —y hg 
for the light quarks is very small and dominated by the charm quark. This 
implies that any interference between the SM and the CMDM is negligible 
except for charm. The new physics effects are thus purely quadratic in both 
the CEDM and the CMDM of the light quarks in our present study. This 
picture fails beyond LO, where qg —y qh proceeds through a heavy quark 
loop 35( , resulting in a more sizable contribution from light quarks which can 
interfere with the new physics. 



The new physics contribution could also be studied in the context of Higgs 
plus one jet (or more) production. This path would have several advantages: 
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1. the NP would be compared to a smaller SM cross-section 

2. the interference between the SM and the CMDM would be included, 
resulting in a linear dependence on NP. 

3. the NP/SM ratio increases with pt so requiring a jet is favorable for the 
purpose of constraining NP. 

At the same time, Higgs plus jets production is a much more complicated 



process [35|, |36| and its study is beyond the scope of this paper. We will 
content ourselves for now with the simpler pp — > hX study, keeping in mind 
that NLO processes may become much more relevant if there are hints for NP. 

In view of these considerations we turn again to MadGraph5 to estimate the 
cross-section pp — > hX due to new physics for different values of the anomalous 
couplings. We begin with d uG at 8 TeV for which we find 

a(pp -> hX) « 2.3 + 0.176 [Re(d uG ) 2 + lm(d uG ) 2 ] pb. (17) 

The SM contribution here, 2.3 pb, arises primarily from cc — > Hg and doesn't 
interfere with the NP which affects only the up-quark. As argued above, we ignore 
this contribution and compare the NP directly with the dominant gluon-fusion pp — > 



hX cross- section for ran = 125.5 GeV [37 



a(pp -> hX) = (19.37 ± 14.7) pb. (18) 

Requiring that the new physics be smaller than the theoretical error leads to the 
constraint 

K|,K| <3.5x lO^GeV -1 . (19) 

We repeat the exercise for 14 TeV, this time including all the light quarks. We 
find 

a(pp -► hX) « 5.58 + 0.7 [Re(d uG f + lm(d uG ) 2 ] 

+ 0.4 [Re(d dG ) 2 + lm{d dG ) 2 ] + 0.1 [Re{d sG f + lm(d sG ) 2 ] 
- 0.034 Re(d cG ) + 0.07 [Re(d cG ) 2 + Im(4c) 2 ] (20) 

This is then compared to the NLO theoretical (dFG) cross-section and error as given 



in Ref. |30j for a mass of M H = 125 GeV 

a(pp -> hX) = (49.97+^) pb (21) 

by requiring the new physics contributions to fall within the la uncertainty. 

The resulting allowed regions, within the la range of Eq. |2TJ and allowing only 
one d q o to be nonzero at a time are shown in Figure [71 For the up, down and strange 
quarks the allowed regions are just circles of larger radii starting with the inner most 
for up ending with the one for strange. The outermost region is for charm and shows 
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FIG. 7: pp — > HX production cross-section at the LHC at 14 TeV as a function of the 
light quarks CEDM and CMDM couplings allowing only one d q c to be nonzero at a time. 
The concentric regions, starting with the smallest one, show the allowed couplings at la 
for u, d, s and c quarks respectively. 



a slight asymmetry due to the interference with the SM. In terms of CEDM and 
CMDM couplings these limits correspond to 



|<|, K| < 1.2 x 1(T 4 GeV" 1 \a% \d 9 d \ < 1.6 x 10~ 4 GeV" 1 

\a 9 s \, \d 9 s \ < 3.3 x 1(T 4 GeV" 1 \d 9 c \ < 3.9 x 10" 4 GeV" 1 

-3.8 x 10~ 4 GeV" 1 < \a a c \ < 4.0 x 10~ 4 GeV -1 



(22) 



We are only aware of existing low energy constraints for the CEDM, where 
the strongest bound arises from the neutron electric dipole moment (dj^). We can 
estimate the contribution of an u p-q uark (down-quark) CEDM to the neutron edm 
using naive dimensional analysis |20l . 138 



d 1 
a N 



as 
—d 9 



(23) 



and from the experimental limit on the neutron edm, d? N < 2.9 x 10 e-cm [21J we 
then find 



d 9 ud < 1.8 x 10" 11 GeV" 1 . 



Similarly from the limit on the A edm, dl < 1.5 x 10 16 e-cm 39j we find 



d 9 < 0.1 GeV" 1 . 



(24) 



(25) 



Using the contribution from the Weinberg operator to the neutron edm, and with 
m c ~ 1.28 GeV, we also obtain 



m c d 9 c < 6 x 10~ 10 . 



(26) 
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V. CONCLUSIONS 

We have studied the constraints that can be obtained on CEDM and CMDM 
type couplings of quarks by studying Higgs boson production at the LHC. For the 
case of the top quark we have compared them to what can be obtained from studying 
top-quark pair production and decay. We find that the pp — >■ tth cross-section can 
provide a significant improvement on the CEDM bound over that from the pp — > ti 
cross-section, and comparable to what can be obtained using T-odd asymmetries 
in the latter. A further improvement is in principle possible by studying T-odd 
correlations in pp — > tth. The low statistics for this channel imply that thousands of 
fb _1 would be needed. The improvement on the top CMDM constraint is modest, 
on the other hand. 

For the other quarks, the Higgs production process may be the only opportunity 
to constrain them at LHC. Our results are summarized in Table [I] In this Table we 

TABLE I: Summary of results for 1<t bounds that can be placed on the CEDM and CMDM 
couplings of quarks at the LHC. The last column shows the effective new physics scale 
than can be probed by the LHC with the given process, the two numbers corresponding 
to the CMDM and the CEDM respectively. 



Process 


CMDM 


CEDM 


A (TeV) 


a(pp ->• tt) 8 TeV 


-0.034 <, m t af <> 0.031 


\m t d 9 t \ <, 0.12 


(1.5, .7) 


a{pp ->■ tt) 14 TeV 


-0.029 <, m t a 9 t <, 0.024 


\m t d?\ ^0.1 


(1.5, .7) 


Ai{pp^tt) 14 TeV 


- 


\m t d 9 t \ ^ 0.009 


(-, 2.5) 


a(pp ->• tth) 14 TeV 


-0.016 <, m t af <> 0.008 


\m t d 9 \ <, 0.02 


(2, 1.7) 


Ai j2 {pp^tih) 14 TeV 


- 


\m t d 9 t \ ^ 0.007 


(-,3) 


a(pp ->• bbh) 14 TeV 


-1.3 x 10~ 4 £ m h a 9 b & 2.4 x 10~ 4 


|m 6 df | ^ 1.7 x 10" 4 


2.7 


a{pp -> hX) 8 TeV 


\al\ <> 3.5 x 10~ 4 GeV" 1 


|<g| ^ 3.5 x 10~ 4 GeV" 1 


1 


a{pp -> hX) 14 TeV 


\a 9 u \ <> 1.2 x 10~ 4 GeV" 1 


\di\ <> 1.2 x 10" 4 GeV" 1 


1.7 


a{pp -> hX) 14 TeV 


\a 9 d \ <> 1.6 x 10~ 4 GeV" 1 


141 ^ 1.6 x 10^ 4 GeV" 1 


1.5 


a(pp -)■ hX) 14 TeV 


|af | <> 3.3 x 10~ 4 GeV" 1 


\d 9 s \ <> 3.3 x 10" 4 GeV" 1 


1 


a(pp -»■ hX) 14 TeV 


\a 9 c \ & 3.9 x 10- 4 GeV" 1 


\d 9 c \ ^ 3.9 x 10- 4 GeV- 1 


1 



show the constraints we obtained on the anomalous color electric and magnetic 
dipole moment couplings of quarks in processes with Higgs bosons (and without for 
the case of top). In the last column we use Eq. [3] to interpret these limits as the 
scale of new physics that can be probed by these processes at LHC under the usual 
assumption that the effective Lagrangian dimensionless couplings d q a are of order 
one. 
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Appendix A: pp — > tt events 

We first use MadGraph5 to compute the ti production cross-section at the LHC at 
8 TeV for a set of values of Im(<i iG ) or Re(d t c) with A = 3 TeV. The results of these 
numerical calculations, including the error estimated by MadGraph5, are plotted in 
Figure [HJ In the same figure we show the result of a quartic fit to these points that 
is guided by theoretical prejudice: we know that the cross-section will be a quartic 
polynomial in the anomalous couplings because they occur once (qq annihilation and 
s-channel gluon fusion diagrams) or twice (t-and-u-chanel gluon fusion diagrams) at 
the amplitude level. We also know that the CP violating coupling Im(o?tG) can 
only contribute to the cross-section through even powers. Finally, from an analytic 
calculation of the squared matrix element for the gluon fusion process, we know 
that the quartic terms in Re(dtc) and in Im(dtG) are the same. All these aspects 
are checked to be satisfied within the statistical error of the event generation (100k 
events per point). The quartic polynomial that fits these points is 

(pp -> tt) © 8TeV 






dtG 



FIG. 8: Cross-section for the process pp — > ti at vS = 8 TeV as calculated with MadGraph5 
for different values of the anomalous coupling dtG and the corresponding quartic fit de- 
scribed in the text. 

1 + 5.37 x 10^ 4 lm(d tG ) 2 + 2.41 x 10" 8 Im(d tG ) 4 + 3.56 x 10" 2 Re{d tG ) 



CSAI 



+ 7.75 x 10" 4 Re(d tG ) 2 + 5.87 x 10" 6 Re(d tG ) 3 + 2.45 x 10~ 8 Re(d tG ) 



(Al) 



We next use MadGraph5 to generate event samples for the process pp — > tt at 
14 TeV for different values of the top CEDM and CMDM. The results of these 



14 



numerical calculations, including the error estimated by MadGraph5, are plotted in 
Figure [9j In the same figure we show a quartic fit to these results, as was described 
above for the 8 TeV case. The resulting fit is given by 



b 



(pp->- tt) @ 14TeV 









/ 








.— . MC 

- Re(d, G ) 
Im(dtc) 


/ 










Jy 











FIG. 9: Cross-section for the process pp — > tt at 14 TeV as calculated with MadGraph5 for 
different values of the anomalous coupling d tG and the corresponding fit. 



a 



&SM 



1 + 6.06 x 1(T 4 lm(d tG ) 2 + 3.28 x 10" 8 lm(d tG ) 4 + 3.58 x 10~ 2 Re(d tG ) 



+ 8.46 x 1(T 4 Re(d tG ) 2 + 6.76 x 10" 6 Re(d tG ) 3 + 3.26 x 10" 8 Re{d tG ) 



(A2) 



Next we generate larger samples of 10 6 events in the di-muon channel, pp — > 
tt — > bbn + pTuu, for 4 values of Im(d tG ) with the cuts of Eq. [8] and the missing $ T > 
30 GeV requirement. These larger samples are necessary to estimate asymmetries 
at the 10 -3 level. We use the four points to fit a linear form for the asymmetry 
since it can only be generated by the interference between CP conserving and CP 
violating amplitudes. The results of our event simulations and the corresponding fit 
are shown in Figure [3j In principle the results could have terms cubic in the CEDM 
but our fits indicate their effect is negligible in this range of Im(d tG ). One would 
expect the asymmetry to be generated by terms of the form Im(dt G ) Re(dt G ) which 
we have not simulated in our MC studies because they are expected at a much lower 
level. We have also checked that several other asymmetries that should vanish (for 
example those linear in the beam momentum) are consistent with zero within the 
statistical error of our event simulation. Finally, we have also verified numerically 
that there are no asymmetries induced by Re(d tG ), as expected. 

In Ref. [14J it was found that (Eq. 13 of that reference but with the notation 
of this paper) Ai rs —4 x 10 -3 \m(d tG ) which is in reasonable agreement with our 
present result but not identical to it. There are two main differences between the two 
calculations: the different set of parton distribution functions; and the MadGraph5 
implementation. Here we implement the Lagrangian into MadGraph5 with the aid 
of FeynRules, whereas in Ref. 14| we directly modified the MadGraph4 code to 
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insert the analytical result for the interference between the SM and the CEDM from 
Ref. 10] . The latter explicitly removes all terms that are not linear in the CEDM. 



Appendix B: pp — > tth events 



We first consider the production of a Higgs boson associated with a top pair, 
pp — > tth, and generate MC samples of 10 5 events at vS = 14 TeV for several values 
of the top CEDM and the top CMDM. We show the results for the cross-sections in 
Figure [TQ1 We show in the same figure the result of a quartic fit to these points for 



(pp^rtth) 8 14TeV 




FIG. 10: Cross-section for the process pp — > tth as calculated with MadGraph5 for different 
values of the anomalous coupling dt G and the corresponding fit. 

extrapolation purposes. The number of parameters in the fit is reduced by noting 
that the cross-section only has quadratic and quartic terms for the CEDM resulting 
in 



a 



&SM 



1 + 1.53 x 10~ 2 lm(d tG ) 2 + 3.32 x 10" 6 lm(d tG ) A + 1.14 x 10" 1 Re(d tG ) 



1.70 x 10" 2 Re(d tG ] 



2.90 x 10" 4 Re(d tG ) 3 + 3.32 x 10" 6 Re(d tG ) 



(Bl) 

The cross-sections obtained with this approximate formula are in agreement with 
the results in Ref. 24J for the linear and quadratic terms. 

We next consider several asymmetries by generating samples of 10 6 events (to 
reach a sensitivity of 10~ 3 ) for di-muon decays of the top pair and using the same cuts 
as in th previous section. The largest asymmetries we find are shown in Figure [51 
In addition we tried the operator 

Oi = p h ■ (p„+ x #,-) (B2) 

but it results in an symmetry too small to measure with 10 6 events. We also tested 
correlations that are known to be zero, such as Pbeam • ij>^+ x P^-), to check that 
we indeed get zero within our statistical uncertainty. In the same manner we found 
that the CMDM does not induce any of the T-odd correlations, as expected. 
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Appendix C: pp — > bbh 

We repeat the exercise of the previous section, this time for the process pp — > bbh. 
To this end we implement the 6-quark CEDM and CMDM into MadGraph5 and 
compute the resulting cross-section for several values of these couplings. The results 
of these runs along with a quartic fit to them is shown in Figure [TTJ The quartic fit 

(pp -> bbh) Q 14TeV 




FIG. 11: Cross-section for the process pp — > bbh in the di-muon channel as calculated with 
MadGraph5 for different values of the anomalous coupling d bG and the corresponding fit. 

to these points enforcing the condition of no odd powers of the CEDM is given by, 



a 



&SRI 



1 + 1.49 x KT 1 Ihi(4g) 2 + 1-02 x lCT 5 Ihi(4g) 4 - 9.76 x 1(T 2 Re(d bG ) 



1.49 x KT 1 Re(d bG ) 2 + 9.73 x 1(T 5 Re(4c) 3 + 1-02 x 10" 5 Re(d bG ) 



Appendix D: pp — > hX 



(CI) 



To constrain the anomalous couplings a 9 and d 9 of the light quarks we consider 
their contribution to inclusive Higgs production. The dominant contribution is from 
the up-quark through the parton level diagrams ug —$■ uh and uu — > gh. Since we 
are not requiring a jet in the final state we have removed the (pt)™™ cut associated 
with the final state quark (gluon) from the MadGraph5 default cuts. We compute 
the resulting cross-section for several values of the couplings along with a quadratic 
fit (because in this case the new coupling appears only once in the amplitudes) for 
the LHC at VS = 8 TeV. 

The quadratic fit we obtain for these points is consistent with having no linear 
term for either the CEDM (as expected from CP) or the CMDM. The latter is 
consistent with the SM contribution being very small in the diagrams with up- 
quarks as it is proportional to the quark mass. The leading SM contribution to the 
tree level processes we evaluate arises from the parton level process cc — > hg and 
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therefore only interferes with the charm quark CMDM. The result of the fit is given 
in Eq. El 

The corresponding results for the LHC at V5 = 14 TeV are shown in Figure [ 
The resulting fit is given in Eq. 







(pp -> hX) 14TeV 


3.5 














i—i MC 




3 




- Re(d» G ) 


/• 
















E 2 
b 




b 1.5 




1 




0.5 


K A 



(pp -ttX) 8 14TeV 




(pp -> hX) 8 14TeV 



(pp -> hX) (Qj 14TeV 





FIG. 12: Cross-section for the process pp — > HX as calculated with MadGraph5 for different 
values of the anomalous coupling d q c and the corresponding fit at y/S = 14 TeV. 
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